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DYNAMICS  OF  AN  INTENSE  RELATIVISTIC  ELECTRON 
BEAM  INJECTED  INTO  FULL  DENSITY  AIR 


I.  Introduction 

When  an  intense  relativistic  electron  beam  (REB)  is  injected  into  the  atmosphere  it  ionizes  the 
molecules  both  in  its  path  and  in  the  surrounding  atmosphere.1  2  The  electric  field  EArj)  produced  by 
the  beam  drives  a  ’‘return"  current  in  this  conducting  plasma  known  as  the  plasma  current,  and  the 
interaction  of  the  injected  electron  beam  with  the  plasma  current  plays  a  major  role  in  determining  the 
stability  and  subsequent  dynamics  of  the  REB  as  it  propagates  through  the  atmosphere.3  Thus  it  is 
important  to  be  able  to  calculate  the  plasma  current  induced  by  a  given  REB  and  to  understand  the 
interaction  between  these  two  currents. 

We  describe  an  experiment  in  which  an  intense  REB  (—1  MV,  —16  kA,  —25  ns)  has  been 
injected  into  the  atmosphere.4  The  induced  plasma  current  reached  a  peak  value  of  almost  1 1  kA  and 
the  motion  of  the  REB  was  studied  using  a  high  speed  streak  camera.  The  intense  REB  was  produced 
using  a  cold  cathode  vacuum  diode  and,  as  in  our  previous  experiments,5  6  was  injected  into  the  atmo¬ 
sphere  through  the  anode  foil. 

The  ionization  of  the  atmosphere  by  the  intense  REB  and  the  generation  of  the  plasma  current 
has  again  been  modeled  using  the  CHMAIR  code7  and  a  simple  circuit  (or  transmission  line)  equation 
in  which  the  REB  is  replaced  by  a  current  source.  The  net  current  predicted  by  this  model  is  compared 
with  the  measured  net  current. 

Manuscript  submitted  on  July  2.  1981. 
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II.  The  Pulse  Generator  and  Diode  System 

The  same  modified  Physics  International  Pulserad  310  pulse  generator  and  cold-cathode,  vacuum 
diode,  as  used  in  the  experiments  described  in  Ref.  5,  were  used  in  the  experiments  presented  here. 
With  the  diode  slightly  mismatched  to  the  40  ft  Blumlein,  the  system  produced  a  single  pulse  of  elec¬ 
trons  inside  the  diode  with  peak  voltage  —  1  MV  and  peak  current  —20  kA:  the  charging  voltage  on  the 
Marx  generator  was  90  kV.  (These  numbers  must  be  compared  to  the  diode  voltage  of  —0.95  MV  and 
the  diode  current  of  —24  kA  in  Ref.  5.) 

The  diode  configuration  is  shown  in  Fig.  1.  It  is  the  same  as  used  in  Ref.  5  except  that  the 
anode-cathode  {A-K)  gap  has  been  increased  from  17  mm  to  25  mm.  The  polycarbonate  prepulse 
switch  was  used  in  series  with  the  carbon  cathode  to  prevent  breakdown  of  the  diode  too  early  in  the 
voltage  pulse.  The  copper  shield  prevented  u-v  preionization  of  the  A-K  gap  and  inhibited  plasma 
formed  at  the  switch  surface  from  expanding  into  the  gap.  The  anode  of  the  diode  was  an  annealed 
titanium  foil,  0.143  mm  thick,  which  bowed  into  the  diode  approximately  5  mm.  The  foil  diameter  was 
10.4  cm  and  the  diode  operating  pressure  was  —  3  x  10“ 5  Torr.  The  diode  current,  IDy  was  monitored 
by  a  A  loop  in  the  wall  of  the  diode  and  the  diode  voltage,  VD,  was  measured  by  a  resistive  divider  at 
the  cathode  shank  and  was  corrected  for  the  voltage  drop  along  the  shank, 

III.  REB  Propagation  in  the  Atmosphere 
Ilia.  Injected  Electron  Beam  Characteristics 

The  characteristics  of  the  electron  beam  injected  into  the  atmosphere  were  measured  using  the 
configuration  shown  in  Fig.  2.  In  this  set  up  the  REB  propagated  inside  a  conducting  tube  made  of  fine 
mesh  (—0.125”)  copper  screen  which  was  attached  to  the  anode  foil  stretcher  at  one  end  and  to  the  foil 
stretcher  of  the  vacuum  Faraday  cup  (VFC)  at  the  other.  The  foil  on  the  VFC  was  titanium  of  the 
same  thickness  as  the  anode  foil  and  the  VFC  was  evacuated  to  a  pressure  of  — 10  m  Torr.  (The  resis¬ 
tive  element  of  this  VFC  is  a  titanium  foil,  — 12  Aim  thick,  which  has  less  inductance  than  the  stainless 
steel  foil  used  previously5  6  and  permits  a  response  time,  rnse  <1  ns.l  The  carbon  collector  of  the  VFC 
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has  a  diameter  of  3.5”  (8.9  cm).  The  cup  was  always  used  with  an  aperture  with  diameter  less  than  3,5” 
to  prevent  the  beam  from  shorting  out  the  insulator  around  the  carbon  collecting  block.  A  fast 
Rogowski  coil  (RC),  /rise  ns,  was  used  simultaneously  inside  the  conducting  tube  to  monitor  the  net 
current  in  the  region  of  beam  induced  ionization.  In  this  configuration  various  graphitic  carbon  aper¬ 
tures  were  placed  in  front  of  the  VFC  to  allow  measurement  of  the  beam  current  distribution,  Ib(r,t), 
from  which  the  beam  current  density  distribution  Jb(r,t)  could  be  obtained.  [The  diagram  shown  in 
Fig.  2  is  not  to  scale— the  actual  distance  between  the  anode  foil  and  the  VFC  was  approximately  6  cm.] 

Figure  3  shows  typical  records  of  the  diode  current,  the  net  current  monitored  in  front  of  the 
aperture  (RCI),  the  net  current  monitored  behind  the  aperture  (RC2),  and  the  beam  current  (meas¬ 
ured  with  the  VFC).  The  diode  current  shows  a  tail  of  current  at  low  voltage  due  to  plasma  formation 
inside  the  diode.  This  "tail"  has  been  observed  previously  in  similar  diode  configurations  (Refs.  5,6) 
and  does  not  correlate  with  the  injected  beam  current:  the  tail  is  not  present  in  /RC1,  /RC2  or  7Vfc- 
Although  /RC2  is  very  similar  in  amplitude  to  Iv fc»  peak  of  7RC2  (and  of  /RC i)  is  delayed  by  —5  ns 
from  that  of  7V fo  showing  that  even  with  the  REB  propagating  inside  a  conducting  tube  there  was 
some  plasma  return  current.  Previously6  we  had  found  that  for  a  beam  with  VD  —  2  MV,  lb  —40  kA, 
and  duration  —50  ns  (FWHM)  propagating  inside  a  similar  conducting  tube  there  was  virtually  no 
plasma  return  current  (/nct  —  lb).  But  in  those  experiments  no  measurement  was  made  of  the  delay 
between  net  and  beam  currents. 

The  measured  result  that  /RC1  >  /RC2,  shows  that  current  was  flowing  at  radii  larger  than  that 
intercepted  by  either  RC2  or  the  VFC.  Therefore  no  direct  measurement  was  possible  of  the  total 
beam  current.  To  eliminate  shot  to  shot  variations  from  the  VFC  measurements  with  different  size 
apertures,  the  VFC  signals  were  normalized  to  /RC1  on  each  shot.  Values  of  /VFc(peak)//RC|  (peak) 
plotted  as  a  function  of  aperture  radius  are  shown  in  Fig.  4.  Notice,  the  shape  of  7Vfc  did  not  change 
with  the  size  of  the  carbon  aperture  as  it  did  in  the  earlier  experiments.6  To  determine  the  total  beam 
current,  the  VFC  measurements  were  extrapolated  to  the  full  anode  foil  radius  of  5.2  cm.  Because  of 
the  uncertainty  in  the  extrapolation  and  the  possibility  of  beam  expansion  between  the  anode  foil  and 
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the  VFC,  the  peak  beam  current  was  determined  as 

Ib  (peak)  -  16  ±  1  kA , 

compared  with  the  peak  net  current  measured  on  RC1  of 

/rci  (peak)  -  14  ±  1  kA. 

Our  best  estimate  of  the  time  dependence  of  the  total  beam  current  was  taken  as  being  the  same  as  that 
of  the  VFC  signal  measured  using  the  largest  aperture,  i.e.,  radius  —3.75  cm  which  is  shown  in  the  bot¬ 
tom  trace  in  Fig.  3. 

Taking  the  VFC  signals  as  Ib(r,t)  the  beam  current  density,  Jb (/■,/),  can  be  determined  for  all  r 
and  t.  The  beam  current  density  corresponding  to  the  time  of  the  peak  total  beam  current  is  shown  in 
Fig.  5.  The  maximum  current  density  was  —0.5  kA/cm2  and  the  radius  containing  —50%  of  the  total 
beam  current  was  —3  cm. 

Clearly  when  studying  the  freely  propagating  electron  beam  no  direct  measurements  could  be 
made  of  the  injected  beam  current  and  current  density.  Therefore  we  kept  the  diode  configuration  and 
generator  parameters  constant  and  assumed  that  the  beam  current  injected  into  the  atmosphere  without 
the  conducting  screen  was  the  same  as  measured  above  with  the  conducting  screen.  As  a  check  on  this 
assumption  we  monitored  the  diode  characteristics,  ID  and  VD ,  on  almost  every  shot.  For  the  above 
measurements 

<ID>  -  20  ±  1  kA 
and 

<  VD>  -  1  ±  0.05  MV, 

whereas  for  the  measurements  on  freely  propagating  electron  beams 

<ID>  -  18  ±  1  kA 
and 

<  VD>  -  1.1  ±  .05  MV. 

The  shapes  of  Id  and  VD  did  not  change  however. 
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Ulb.  Characteristics  of  the  Freely  Propagating  Electron  Beam 

In  the  free  beam  configuration  shown  in  Fig.  6,  the  Rogowski  coil,  shielded  from  beam  electrons 
by  the  foil  stretcher,  monitored  the  net  current,  lnel(t)y  on  all  shots.  To  eliminate  shot  to  shot  varia¬ 
tion  the  average  net  current  was  obtained  as 

I  (peak) 

</nel(peak)>  -  <  - rr->  '  </D(peak)>  -  7  kA. 

/fltpeaK; 

Typical  oscilloscope  traces  for  //>,  Fd,  and  /net  are  shown  in  Fig.  7.  Clearly  the  net  current  profile  is 
much  broader  than  that  of  lD  and  the  peak  of  /net  is  delayed  from  the  peak  of  ID,  both  of  which  indi¬ 
cate  significant  plasma  return  current. 

The  time  integrated  photograph  shown  in  Fig.  6  suggests  that  the  freely  propagating  electron  beam 
disrupts  after  a  region  (—50  cm)  of  quasi-stable  propagation.  In  this  region  coherent  betatron  oscilla¬ 
tions  can  be  seen  as  they  were  in  our  earlier  high  current  experiment.6 

IIIc.  Diode  Asymmetry  Effects 

Given  the  symmetric  position  of  the  electron  beam  in  the  screened  enclosure  shown  in  Fig.  6,  we 
fully  expected  that  the  disruption  of  the  beam,  after  its  initial  period  of  quasi-stable  propagation,  would 
be  random.  Therefore  we  were  surprised  when  in  a  sequence  of  approximately  twenty  consecutive 
shots  the  electron  beam  consistently  deflected  downwards  as  it  disrupted.  Subsequently  we  discovered 
an  asymmetry  in  the  diode  and  found  that  it  was  responsible  for  the  non-random  disruption  of  the 
beam.  This  asymmetry  is  illustrated  in  Fig.  8.  The  copper  shield  had  inadvertently  been  placed  with  its 
axis  of  symmetry  a  full  2  mm  below  the  axis  of  the  cathode.  The  asymmetry  was  reversed  and  the 
beam  deflected  upwards.  Typical  results  caused  by  this  asymmetry  are  shown  in  Fig.  9. 

The  faulty  copper  shield  was  replaced  with  one  that  centered  on  the  cathode  assembly  so  that  the 
whole  assembly  was  coaxial  within  —  ±0.2  mm.  Thereafter  the  electron  beam  frequently  propagated 
further  (—1  m)  before  disrupting  and  deflected  randomly. 
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Effects  that  are  believed  to  have  been  caused  by  asymmetries  inside  diodes  have  apparently  been 
seen  before8  but  no  detailed  information  was  available.  While  questioning  how  asymmetries  inside  a 
diode  might  affect  beam  propagation  outside  the  diode,  we  realized  that  with  an  extended  cathode  such 
as  the  one  we  have  used,  electron  emission  may  not  be  uniform.  Then  if  the  point  of  emission  on  the 
cathode  were  to  start  at  a  given  position  in  relation  to  the  asymmetry  and  to  move  during  the  beam 
pulse  to  a  position  diametrically  opposite  that  initial  position,  the  later  parts  of  the  beam  pulse  would 

mostly  be  emitted  from  the  anode  foil  displaced  to  one  side  of  the  earlier  parts  of  the  beam  pulse. 

Then  because  of  the  time  required  for  atmospheric  conductivity  to  be  established  and  to  decay,  the 
plasma  return  currents  would  tend  to  be  displaced  with  respect  to  the  beam  current.  But  this  is  an 
unstable  situation  in  which  any  initial  small  displacement  of  the  beam  is  amplified  as  it  convects  back 
into  the  beam  pulse. 

Hid.  Time  History  of  Fluorescence  Radiation 

In  preparation  for  the  intended  study  of  the  dynamics  of  the  intense  REB  propagating  in  the 

atmosphere  using  fast  streak  photography,  a  careful  measurement  was  made  of  the  time  history  of  the 

induced  fluorescence  radiation  and  its  correlation  with  the  beam  current.  The  streak  camera  (an  Ima- 
con  790)  was  known  to  contain  an  S20  photocathode  on  a  1/4  inch  glass  window;  therefore,  to  measure 
the  time  history  of  the  fluorescence  radiation,  we  used  a  IP28  photomultiplier  and  a  glass  lens.  Thus 
both  the  streak  camera  and  the  photomultiplier  responded  to  approximately  the  same  wavelength  range, 
X  ^3500  A.  The  spectrum  of  fluorescence  radiation  from  a  comparable  REB  interacting  with  the  atmo¬ 
sphere  had  already  been  observed6:  it  consists  primarily  of  emission  in  the  second  positive  system  of 
;V2.  In  the  range  X  ^3500  A  most  of  the  emitted  radiation  still  comes  from  the  N2  2P  system  though 
the  strongest  emission  in  that  system  comes  at  shorter  wavelengths  [N2  2P  (0-0)  at  X  3371  A]. 

Because  measurements  of  the  time  history  of  the  fluorescence  radiation  were  made  on  the  freely 
propagating  REB,  comparison  could  not  be  made  directly  with  the  injected  beam  current,  i.e.  /Vfo 
instead  comparison  was  made  with  the  diode  current,  ID.  Earlier  measurements  (Fig.  3)  have  already 
shown  that  /VFc  and  b  rise  simultaneously  within  ±  1  ns.  Measurements  were  made  viewing  a  5  cm 
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length  of  the  beam  path,  looking  perpendicularly  to  that  path  at  positions  7  cm  and  50  cm  from  the 
anode  foil.  The  photomultiplier  cable  length  was  adjusted  so  that  a  light  signal  emitted  at  the  anode  foil 
was  observed  on  the  oscilloscope  trace  as  simultaneous  within  ±  1  ns  with  the  signal  on  the  ID  cable. 
The  photomultiplier  was  completely  enclosed  in  a  1  inch  thick  lead  housing  to  stop  the  scattered  x- 
radiation  (sky  shine)  with  an  additional  2  inch  thick  lead  wall  between  the  photomultiplier  and  the 
diode  to  stop  the  direct  x-radiation.  Light  input  to  the  photomultiplier  was  reflected  from  a  small  mir¬ 
ror  inside  the  lead  housing. 

With  the  entrance  to  the  photomultiplier  blocked  and  the  lead  wall  removed,  energetic  x-rays, 
created  mostly  in  the  anode  foil,  reached  the  photomultiplier.  This  x-ray  signal  was  of  shorter  duration 
(FWHM  —20  ns)  than  ID  (FWHM  —30  ns),  had  a  faster  rise-time  (—5  ns)  and  was  delayed  from  lD 
by  approximately  5  ns.  Thus  the  x-ray  signal  corresponded  exactly  to  that  part  of  the  diode  current  for 
which  the  diode  voltage  was  ^0.5  MV,  which  is  very  reasonable  since  the  photomultiplier  was  still 
enclosed  in  the  1  inch  thick  lead  housing.  [Note  that  because  of  the  different  path  lengths  travelled  by 
the  optical  and  x-radiation,  the  x-ray  signal  on  the  photomultiplier  was  actually  recorded  on  the  oscillo¬ 
scope  —  10  ns  before  lD .]  With  the  lead  wall  in  place  no  appreciable  x-ray  signal  was  observed. 

A  typical  photomultiplier  signal  showing  the  fluorescence  radiation  is  shown  in  Fig.  10.  The  pho¬ 
tomultiplier  signal  has  a  slower  rise-time  (—15  ns)  than  1D  (—10  ns),  is  delayed  from  ID  by  approxi¬ 
mately  6  ns,  but  has  approximately  the  same  width  as  ID.  [Actually  the  photomultiplier  signal  shows 
some  distortion  due  to  electrical  noise  and  is  probably  a  little  shorter  than  lD  as  is  the  beam  current 
(FWHM  —27  ns).]  The  photomultiplier  signals  at  the  two  different  axial  positions  were  the  same 
within  experimental  reproducibility. 

For  these  measurements  we  have  used  the  same  photomultiplier  system  as  we  used  in  the  earlier 
measurement6  when  we  found  that  the  light  output  was  synchronous  with  the  diode  current  within  ±5 
ns.  However  for  these  measurements  we  have  taken  more  care  to  match  cable  and  light  path  lengths, 
and  checked  this  matching  using  an  ultra  fast  light  source.  The  rise-time  of  the  photomultiplier  system 


was  <3  ns  which,  together  with  the  electrical  noise  on  the  signals,  limited  the  measurement  accuracy 
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to  approximately  ±  1  ns  and  would  have  prevented  detection  of  the  difference  in  the  time  of  emission 
for  the  two  axial  positions  even  for  signals  with  shorter  rise-times.  On  the  other  hand  the  extended 
rise-time  of  the  light  signal  compared  to  that  of  the  diode  current  and  the  delay  between  the  current 
signal  and  the  light  signal  were  both  readily  measured. 

Since  most  of  the  radiation  seen  by  the  photomultiplier  (and  by  the  streak  camera)  comes  from 
the  second  positive  system  of  N2 ,  the  intensity  observed  is  directly  proportional  to  the  population  of  the 
upper  state,  the  N2(C)  state: 

p 0-1  a  ^0-1  N 2^0 

where  Pq_ t  is  the  power  radiated  per  unit  volume  in  the  2P(0—  1)  band,  is  the  Einstein  transition 
probability  for  that  band,  and  N2{C)  represents  the  population  density  of  the  jV2(C)  state.  But 

~  { N 2 ( C ) }  “  t)  -  (£ A  +  qN2)  NAC) 

at  ax  e 

where  7)  is  the  efficiency  for  producing  N2(C)  states  by  beam  excitation,9  — 10~4  eV~l;  dE/dx  is  the 
energy  deposited  per  cm  of  path  length  by  a  beam  electron;10  is  the  total  Einstein  transition  proba¬ 
bility  for  the  N2(C)  state,11  ~  2.5  x  107  sec-1;  and  q  is  the  quenching  rate  for  collisional  depopulation 
of  the  N2(C)  state,11  ~  1.1  xlO-11  cm3  sec-1.  At  ambient  pressure  and  temperature 

qN2  —  2.2  x  108  sec-1; 

therefore  the  equilibrium  concentration  of  the  N2(C)  state  is  given  by 

(0  ^  . 
eqN  2 

The  time  taken  to  reach  this  equilibrium  concentration  is  approximately, 

r  - - 77 - -  10-8  sec. 

qN  2 

which  results  in  a  delay  between  the  switching  on  of  the  electron  beam  and  the  emission  of  radiation  of 
approximately  y r  (or  5  x  10”9  sec)  for  the  case  of  a  beam  pulse  whose  duration  is  longer  than  r. 

Thus  the  measured  delay  of  6  ±  1  ns  is  in  excellent  agreement  with  the  delay  predicted  for  beam 
induced  fluorescence  radiation. 


8 


NRL  MEMORANDUM  REPORT  4557 

Ille.  Streak  Photography  and  Beam  Dynamics 

As  indicated  earlier,  a  major  objective  of  this  experiment  was  to  study  the  dynamics  of  an  intense 
REB  propagating  in  the  atmosphere  using  fast,  time-resolved  photography.  The  camera  used  was  an 
Imacon  model  790,  image  converter  camera  which,  operating  in  its  streak  mode,  allows  writing  speeds 
of  up  to  1  ns/mm  on  the  film  plane.  The  photocathode  on  this  camera  is  only  —20  mm  tall  and  the 
camera  streaks  horizontally.  Therefore  to  observe  the  typical  beam  displacements  indicated  in  the 
open-shutter  photographs  (—  ±50  cm),  the  image  on  the  photocathode  had  to  be  demagnified  approxi¬ 
mately  50:1.  The  shortest  focal  length  lens  available  that  was  capable  of  this  field  of  view,  was  a  small 
achromat  of  aperture  f/5  and  focal  length  50  mm.  This  lens  was  used  for  streak  photography  at  all  axial 
positions  except  ;  =  7  cm  where  an  f/2,  15  cm  focal  length  lens,  which  gave  a  demagnification  of 
—  7:1,  was  used.  To  avoid  electrical  pick-up,  which  caused  pre-firing,  the  streak  camera  was  contained 
in  a  screened  enclosure.  It  was  placed  approximately  1  m  from  the  diode  and  viewed  the  electron  beam 
through  a  series  of  mirrors  (including  one  rear-surfaced  lucite  mirror,  1/4"  thick)  for  all  axial  positions 
except  z  *  7  cm  where  a  direct  view  was  taken.  The  output  of  the  streak  camera  was  recorded  on 
Polaroid  Type  47  film  and  the  streak  speed  was  calibrated  against  a  fast  signal  generator. 

Typical  streak  records  are  shown  in  Figs.  11  and  12  for  axial  positions  of  z  «=  60  cm  and  z  **  88 
cm  respectively.  The  actual  writing  speed  on  these  records  was  1.6  ns/mm  which  because  of  the  width 
of  the  slit  on  the  streak  camera  gave  a  time  resolution  of  —1.2  ns.  [The  length  observed  along  the 
electron  beam  corresponded  to  Az  <  4  cm.]  Thus  on  these  photographs  the  time  resolution  is  actually 
limited  to  —5  ns  by  the  effective  rise-time  of  the  radiation  being  observed  (see  Section  Illd).  Although 
not  readily  visible  in  these  reproductions  (Figs.  11  and  12),  in  the  original  data,  the  characteristic 
trumpet  shape  of  the  head  of  the  REB  was  clearly  seen.  The  width  of  the  trumpet  was,  however,  —4 
ns  and  thus  was  limited  by  the  radiative  relaxation  time. 

From  the  streak  photographs  it  was  clear  that: 

(a)  The  REB  always  remained  as  a  single  current  filament,  with  nearly  constant  diameter  except 
at  the  beam  head. 
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(b)  The  REB  was  emitted  from  the  diode  approximately  on  axis. 

(c)  The  REB  frequently  deviated  from  the  diode  axis  first  in  one  direction  then  in  the  opposite 
direction. 

By  assuming  that  the  REB  propagated  away  from  the  diode  with  a  constant  axial  velocity, 

v2-0c. 

we  were  able  to  interpret  the  streak  photographs  of  a  sequence  of  shots  in  which  the  cathode  asym¬ 
metry  produced  a  nearly  reproducible  beam  deflection  as  a  set  of  snapshots  of  the  electron  beam  over 
the  axial  range  0<  z  ^  108  cm.  This  analysis  is  shown  in  Fig.  13,  where  the  dashed  lines  represent 
that  part  of  the  REB  still  to  be  emitted  from  the  diode  and  the  vertical  line  at  z  —  380  cm  represents  the 
nominal  range  of  1  MeV  electrons  in  the  atmosphere.  Clearly  the  deviations  of  the  REB  from  its  axis 
are  too  large  for  v:  to  be  constant  but  the  data  does  not  merit  any  more  complex  interpretation. 

IV.  Modeling  REB  Propagation  in  the  Atmosphere 

In  our  previous  reports5,6  reference  was  made  to  "circuit  models"  that  were  used  to  calculate  the 
net  current,  /net(/),  for  a  given  injected  beam  current,  40),  but  no  details  of  the  models  were  given. 
In  this  section  we  describe  in  detail  the  two  models  used,  (a)  the  lumped-parameter  circuit  model,  and 
(b)  the  transmission-line  circuit  model,  and  use  the  models  to  predict  7netO)  in  the  situations  described 
in  Section  III.  These  models  are  shown  schematically  in  Figs.  14  and  15,  respectively.  In  both  models 
we  assume  that  the  REB  propagates  within  a  grounded  return  conductor  with  constant  radius,  Rs,  and 
that  the  beam  is  uniform  with  radius  rb.  Then,  as  we  shall  see,  model  (a)  is  relevant  to  the  situation 
shown  in  Fig.  2,  and  Fig.  14a,  where  the  REB  is  stopped  by  a  beam  collector,  the  vacuum  Faraday  Cup, 
which  is  connected  to  the  diode  ound  plane  by  a  coaxial  return  conductor.  Model  (b)  is  relevant  to 
the  case  of  the  REB  propagating  freely  in  the  atmosphere  within  some  much  larger  screened  enclosure 
as  depicted  in  Fig.  6  and  Fig.  15a. 

C  .sider  the  electron  beam  as  it  exits  the  grounded  diode  (located  at  z  —  0)  and  continues  pro¬ 
pagating  within  the  beam-centered,  grounded  return  screen,  /?s(a)  —7.5  cm  and  Rs( b)  —150  cm. 
Assuming  axial  symmetry  is  preserved.  Maxwell's  equations 
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vxfi.iiijj  +  afj  +  iM 
c  c  8r 

Vxf.-iM 

c 


lead  to  an  on-axis  electric  field  given  by 

8 


a2^ 

2 


£.--f(I./„) -L  f 

8z  c2  J°  r  8/ 

where  the  net  current 

/,-  J0*’  <*■  2ir  r  (J^+tr  £:) 

and  the  electrostatic  potential 

r°°  c 

^  “  Jo  rfr  E'  "  J-«  *  Jo  *  2lr 

with  the  charge  density  p  satisfying 


p{r,z',t\ 


V  r2  +  ~Tz‘  —  z)2 


i£. 


+  V  •  (A  +  o-  £)  -  0 


3f  '  '  6 

p(r,  -  z,t )  -  -  p(r,  +  z,t). 

Here,  Jb  is  the  beam  current  density,  or  is  the  air  plasma  conductivity,  and  the  distributed  inductance 

Uiz,+<rE.) 


z.. .  4 /;■*/> 

c2  J o  r  J 0  / 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


The  displacement-current  terms  can  be  dropped  once  the  air  conductivity  satisfies 

47rcr  (/)  •  t  >  ( RJrb )2.  (7) 

This  condition  is  typically  satisfied  within  a  few  nsec  after  the  start  of  the  beam  pulse.  Eq.  (1)  then 

reduces  to 


C---T,  <M.>  -  <*> 

Once  condition  (7)  is  satisfied,  any  excess  charge  within  the  plasma  is  "instantaneously"  tran¬ 
sported  to  the  plasma  edges.  This  charge  generates  electrostatic  racial  fields  which  can  cause  air  break¬ 
down  and  thereby  extend  the  effective  plasma  radius  to 


Rc  >  rb. 

A  crude  estimate  for  R t.  is  afforded  by  the  breakdown  condition 


(9) 
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Rc-  ~^r 
PcEb 


(10) 


where  the  air  breakdown  field  Eb  =  30  kV/cm.  The  electrostatic  potential  <t>  can  then  be  estimated  by 
assuming  that  a  distributed  plasma  charge. 


3  z  Jo 


o  dr2lT 


rp. 


resides  at  radius  Rc  (within  a  thin  sheath  of  thickness  8  «  Rc).  From  Eqs.  (2)  and  (4), 


_3_ 

dr 


1 0. 

bz 


bz 


(11) 


(12) 


Note  that  image  charges  reside  at  Rs. 


The  problem  is  complicated  by  a  mixture  of  boundary  conditions:  (i)  the  beam  current  lb  is 
invariant  in  the  beam  frame,  traveling  at  velocity  j3e,  while  (ii)  the  net  current  satisfies  in  the  lab  frame 


bz 


-  0. 


(13) 


This  mixture  typically  generates  both  a  transient  solution,  which  decays  in  a  time  <  21.//?.,  and  a  non- 
transient  solution,  which  persists  for  a  time 


(14) 


(15) 


h  m  +  LJR:, 

where  r  p  is  the  pulse  duration  and  the  distributed  plasma  conductance 

r*. 

R.  1  —  I  dr  2ir  ra- . 

*  •'0 

Initially,  the  plasma  conductance  is  low,  which  suggests  that  the  transient  solutions  can  often  be 
ignored.  We  thus  focus  only  on  the  non-transient  solutions,  which  generally  persist  much  longer. 
(Note  that  ignoring  the  transients  is  justified  even  if  high  Rz~~l  were  initially  present  due  to  a  preionized 
channel,  providedr p  »  L.J R:.) 

Consider  first  a  configuration  in  which  a  grounded  collector  plate  (e.g.,  a  Faraday  cup)  or  screen  is 
placed  upstream  from  the  diode  a  short  distance 


z0  «  /3Jcr, 


(16) 
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where  r ,  is  the  beam  rise  time.  This  plate  or  screen,  plus  the  grounded  diode  plate  at  z  -  0,  force  the 
radial  electric  fields  to  zero,  and  thereby  effectively  isolate  each  passing  beam  segment  from  long-range 
upstream  space  charge.  The  non-transient  solutions  are  then  inductively  generated  and  depend  only  on 
( rj ),  with 

(17) 

This  latter  equation  represents  the  lumped-parameter  circuit  limit,  for  which  distributed  capacitive 
effects  disappear.  A  similar  limit  arises  if  the  screen  radius  is  reduced  to 

Rs^Rc,  (17) 

or  if  the  screen  and  charge-layer  radii  satisfy 

Rs  >  Rc  »  Per  r.  (18) 

In  the  absence  of  upstream  ground  screens  or  plates,  and  assuming  that 

P ct r  >  Rs  >  Rc,  (19) 

the  nontransient  solutions  can  be  derived  from 

-  -7-  U  /„)  (20) 

dr 

where  r  -  /  -  z/)3cand  the  reduced  distributed  inductance 

L  -  L:  -  1  /P2c2C:  (21) 

with  distributed  capacitance 

C.  -  [ln(/?s2/*c2)]-'.  (22) 

Note  that  the  temporal  profile  of  the  nontransient  net  current  /„  would  still  be  independent  of  z,  except 
for  a  time  shift  corresponding  to  the  position  zof  observation.  It  is  also  interesting  to  note  that  Lz  can 
be  expressed  as  the  sum  of  the  inductance,  I0*  out  to  radius  Rc ,  plus  the  inductance  from  Rc  to  Rs , 

which  equals  -V  (Rs2/Rc2)-  Hence,  in  the  limit  p 2  —  1,  the  effective  inductance 

<r 

L-L0;  (23) 

i.e.,  the  effective  screen  radius  for  the  nontransient  solutions  becomes  equal  to  the  charge-layer  radius 
Rc>  while  the  actual  screen  radius  /?,  >  Rc  becomes  irrelevant. 
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The  equivalent  circuits  representing  these  two  limits,  (a)  the  lumped-parameter  circuit,  and  b)  the 
transmission  line  circuit  are  shown  in  Fig.  14b  and  Fig.  15b,  respectively.  In  our  prevous  reports5  6  we 
used  the  lumped-parameter  circuit  model  to  represent  both  the  REB  propagating  into  a  Faraday  cup  and 
the  freely  propagating  REB.  But  in  this  report  we  shall  use  each  model  in  its  proper  limit. 

Taking  the  beam  parameters  from  Section  III  we  approximate  the  injected  REB  by  a  sine-wave, 

IbU)  -  lb  sin  (< ot)  for  0  <  <  it 

where  lb  *  16  ±  1  kA  and  a >  *  tt/40  x  10~9  s'1.  This  idealized  beam  current  is  shown  in  both  Fig. 
16  and  Fig.  17  as  the  dashed  curve,  and  compares  closely  to  the  measured  beam  current  which  is  taken 
as  / vfc^ )  normalized  to  16  kA.  The  beam  current,  shown  in  Fig.  4,  is  taken  as  a  Bennett  profile  for 
r  ^  3cm  and  as  a  linear  function  of  radius  for  r  >  3  cm.  The  beam  radius  rb  is  essentially  3  cm. 

With  this  input  data  and  using  the  CHMAIR  code7  to  determine  the  conductivity  of  the  air  over 
the  beam  profile,  the  circuit  equations  of  the  two  equivalent  circuits  have  been  solved  numerically  to 
yield  the  net  currents  shown  as  solid  curves  in  Figs.  16  and  17.  In  these  figures  measured  average  net 
currents  (x)  have  been  normalized  to  the  calculated  peak  net  currents.  The  actual  average  values  of  the 
measured  /nct(peak)  are  also  shown. 

As  can  be  seen  the  agreement  between  calculated  and  measured  net  currents  is  really  rather  good, 
especially  if  one  realizes  that  the  real  beam  current  does  not  go  to  zero  at  t  -  40  ns  but  has  a  tail 
approximately  equal  to  the  discrepancy  between  the  two  net  currents.  For  the  freely  propagating  beam 
the  measured  peak  net  current  was 

^  Ior\  ^ 

- -  0.5  ±  0.07 

while  the  value  predicted  from  the  transmission  line  equivalent  circuit  model  was 

/nei  (peak) 

— - 0.6. 

V.  Conclusions 

We  have  described  experiments  in  which  an  intense  relativistic  electron  beam  was  injected  into 
the  normal  atmosphere.  The  beam  characteristics  were  peak  voltage  —  1  MeV  and  peak  current  —16 
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kA,  and  the  time  variation  of  the  beam  current  was  well  represented  by  a  simple  half  sine-wave  with  <o 

—  0.8  x  108  s'1.  The  radial  distribution  of  the  beam  current  was  approximately  a  Bennett  distribution 
with  radius  —3  cm,  except  that  outside  the  Bennett  radius  the  beam  current  varied  linearly  with  radius. 
The  maximum  beam  current  density  was  — 0.55  kA/cm2. 

This  beam  was  observed  as  it  propagated  through  the  atmosphere  in  visible  light  caused  by  direct 
beam  excitation  of  air  molecules.  Using  an  Imacon  790  Image  converter  camera  as  a  streak  camera 
with  all  glass  optics,  the  time  resolution  of  —4  ns  was  limited  by  the  quenching  rate  of  the  N2 
molecules.  Otherwise  the  time  resolution  would  have  been  —  1  ns  which  suggests  that  with  the  more 
sensitive  new  cameras  (e.g.,  the  Imacon  500)  a  time  resolution  of  <0.5  ns  should  be  attainable  using 
radiation  emitted  in  the  first  negative  system  of  the  molecular  ion,  Vf  (Vf  1 . V  ( 0  -  0)  at  wavelength 

—  3914  A].  The  streak  photographs  show  that  the  REB  remained  as  a  single  current  filament  which 
"oscillated”  about  its  initial  direction  as  it  propagated  through  the  atmosphere,  and  that  while  the  "body" 
of  the  beam  was  confined  to  a  radius  of  <5  cm,  the  "nose"  was  expanded  to  give  the  characteristic 
trumpet-like  shape. 

Injection  of  this  relatively  fast-rising  intense  REB  into  the  atmosphere  drove  considerable  return 
current  in  the  beam-produced,  air  plasma.  This  return  current  was  determined  by  measuring  the  net 
current  flowing  in  the  beam  and  was  calculated,  for  the  freely  propagating  REB,  using  a  transmission 
line  circuit  model.  The  peak  net  current  was  approximately  7  kA  but  this  peak  occurred  well  after  the 
peak  beam  current.  The  peak  return  current  was  almost  11  kA  and  occurred  close  to  the  peak  of  the 
beam  current.  The  net  current  was  calculated  from  the  transmission  line  circuit  model  coupled  to  the 
CHMAIR  code  to  calculate  the  beam  induced  air  conductivity.  It  agreed  within  the  experimental  uncer¬ 
tainties  (—  ±  10%)  with  the  measured  net  current. 

While  performing  these  experiments  we  have  found  that  asymmetries  of  as  little  as  —2  mm  in  the 
cathode  assembly  cause  the  emitted  electron  beam  to  deviate  from  its  straight  axial  path  within  a  few 
tens  of  centimeters  from  the  anode  foil.  The  azimuthal  direction  of  this  deviation  is  directly  related  to 
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that  of  the  cathode  asymmetry.  Such  an  effect  was  not  anticipated  especially  for  such  small  asym¬ 
metries.  The  diameter  of  the  carbon  cathode  was  25  mm  and  the  diameter  of  the  copper  shield  behind 
it  was  —57  mm,  while  the  inside  diameter  of  the  coaxial  ground  conductor  of  the  diode  assembly  was 
305  mm.  Although  the  effect  caused  by  this  asymmetry  was  clearly  identified,  no  details  of  the 
mechanism  involved  have  been  determined. 
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Fig.  1  —  The  Diode  configuration 
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DIODE  CURRENT 


PROBE 

Fig.  2  —  Configuration  for  measuring  the  injected  beam  current,  /b,  as  a  function  of  radius  and  time 
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ROGOWSKI  COIL 


CATHODE 


ANODE  FOIL' 


H  K 

10  cm 


/COPPER  SCREEN  ENCLOSURE 


Configuration  for  studying  a  "freely"  propagating  REB  injected  into  a  large  screen  room  (3x3x5m) 
containing  STP  air.  Diode  center  line  is  1.7  m  from  ceiling  and  1.3  m  from  floor 
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Fig.  7  —  Electron  beam  characteristics  for  a  freely  propagating  beam.  Diode 
vot-age  (Vq).  diode  current  f/D),  anticipated  beam  current  Ub)  taken  as  the 
current  previously  measured  using  the  vacuum  Faraday  cup  (/v fc  ^ r  “  5.2)) 
and  the  net  current  /n€t  measured  with  the  Rogowski  coil. 
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Fig.  10  —  Comparison  of  diode  current  and  PM  tube  signal,  /PM,  for  white  light  emission 
produced  by  beam  excitation  of  air  molecules. 


;n  shutter  and  streak  photographs  of  electron  beams  propagating  freely  in  the  atmosphere; 
pictures  (on  right)  show  the  beam  60  cm  from  the  diode,  traveling  to  the  right. 


Fig.  )  )(b)  —  An  enlargement  of  the  streak  photograph  shown  top  right  in  (a), 
but  printed  at  high  contrast  to  bring  out  the  trumpet  shaped  nose. 


Fig  13  —  A  representation  of  the  RFB  propagating  "freely"  in  air  as  deduced  from  streak  photographs 
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RETURN  CONDUCTOR 


a) 


b) 


Fig.  14  —  Equivalent  circuit  for  nontransient  solutions  in  the  limit  of  the  lumped-parameter  circuit  model.  <a)  A  schematic 
diagram  showing  an  REB  with  radius  rb  propagating  in  a  coaxial  return  conductor  with  radius  R$  which  ends  on  a  beam  col¬ 
lector.  (b)  The  equivalent  lumped-parameter  circuit. 


Fig.  15  —  Equivalent  circuit  for  non-transient  solutions  in  the  limn  of  the  transmission  line  circuit  model,  (a)  A  schematic 
diagram  showing  the  REB  with  radius  rb  propagating  in  a  coaxial  return  conductor  of  radius  Rs  which  continues  indefinitely 
The  atmosphere  around  the  REB  becomes  ionized  to  a  radius  at  which  there  is  a  charge  layer,  (b)  The  equivalent 
transmission  line  circuit. 


FIORITO.  FORDHAM.  GREIG.  PECHACEK,  SETHlAN,  FERNSLER.  AND  HALLE 


Oi 

s 

p 


1N3HHHO 


o 


b 


Fig.  16  —  The  net  current  calculated  using  the  lumped  parameter  circuit  model  for  a  REB  propagating  in  full  density  air  inside  a  coaxial 
conducting  lube.  •  /VK  measurements  normalized  to  16  kA,  x  /RCJ  measurements  normalized  to  13  k A,  —  a  sine  wave  used  to  model 
/b>  and - the  calculated  net  current.  4>  shows  the  average  value  of  the  measured  peak  net  current. 


big.  17  _  The  net  current  calculated  using  the  transmission  line  model  for  a  RbB  propagating  freely  in  full  density  air.  •  /VFC  measure¬ 
ments  normalized  to  16  kA,  x  /Rn  measurements  normalized  to  9.6  kA,  —  a  sine  wave  used  to  model  /b,  and - the  calculated  net 

current.  ‘I»  shows  the  average  value  of  the  measured  peak  net  current. 
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